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Abstract 

We study the phase space geometry associated with index 2 saddles of a potential energy surface 
and its influence on reaction dynamics for n degree-of-freedom (DoF) Hamiltonian systems. In 
recent years similar studies have been carried out for index 1 saddles of potential energy surfaces, 
and the phase space geometry associated with classical transition state theory has been elucidated. 
In this case the existence of a normally hyperbolic invariant manifold (NHIM) of saddle stability 
type has been shown, where the NHIM serves as the "anchor" for the construction of dividing 
surfaces having the no-recrossing property and minimal flux. For the index 1 saddle case the stable 
and unstable manifolds of the NHIM are co-dimension one in the energy surface and have the 
structure of spherical cylinders, and thus act as the conduits for reacting trajectories in phase space. 
The situation for index 2 saddles is quite different, and their relevance for reaction dynamics has not 
previously been fully recognized. We show that NHIMs with their stable and unstable manifolds 
still exist, but that these manifolds by themselves lack sufficient dimension to act as barriers in the 
energy surface in order to constrain reactions. Rather, in the index 2 case there are different types 
of invariant manifolds, containing the NHIM and its stable and unstable manifolds, that act as co- 
dimension one barriers in the energy surface. These barriers divide the energy surface in the vicinity 
of the index 2 saddle into regions of qualitatively different trajectories exhibiting a wider variety 
of dynamical behavior than for the case of index 1 saddles. In particular, we can identify a class 
of trajectories, which we refer to as "roaming trajectories", which are not associated with reaction 
along the classical minimum energy path (MEP). We illustrate the significance of our analysis of 
the index 2 saddle for reaction dynamics with two examples. The first involves isomerization on 
a potential energy surface with multiple (four) symmetry equivalent minima; the dynamics in the 
vicinity of the saddle enables a rigorous distinction to be made between stepwise (sequential) and 
concerted ("hilltop crossing") isomerization pathways. The second example involves two potential 
minima connected by two distinct transition states associated with conventional index one saddles, 
and an index two saddle that sits between the two index one saddles. For the case of non-equivalent 
index one saddles, our analysis suggests a rigorous dynamical definition of "non-MEP" or "roaming" 
reactive events. 

PACS numbers: 05.45.-a, 45.10.Na, 82.20.Db, 82.30.Qt 
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I. INTRODUCTION 



Transition state theory has long been and continues to be a cornerstone of the theory 
of chemical reaction rates l|, [2, l^, O, O]. In a number of papers in recent years it has 
been shown that index one saddles [7|] of the potential energy surface give rise to a variety 
of geometrical structures in phase space, enabling_^the realization of Wigner's vision of a 



transition state theory constructed in phase space j8| 
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Following these studies, it is natural to investigate the phase space structure associated 
with saddles of index greater than one, and to elucidate their possible dynamical significance. 
In this note we describe the phase space structures and their influence on transport in phase 
space associated with index two saddles of the potential energy surface for n degree-of- 
freedom (DoF) deterministic, time-independent Hamiltonian systems. 

The phase space manifestation of an index one saddle of the potential energy surface is an 
equilibrium point of the associated Hamilton's equations of saddle-center-. . .-center stability 
type. This means that the matrix associated with the linearization of Hamilton's equations 
about the equilibrium point has one pair of real eigenvalues of equal magnitude, but opposite 
in sign (±A) and n — 1 pairs of purely imaginary eigenvalues, iicufc. A; = 2, . . . , n. 

The phase space manifestation of an index two saddle of the potential energy surface is an 
equilibrium point of the associated Hamilton's equations of saddle-saddle-center-. . .-center 
stability type. The matrix associated with the linearization of Hamilton's equations about 
the equilibrium point then has two pairs of real eigenvalues of equal magnitude, but opposite 



3Q|. 



in sign (±Ai, ±A2) and n — 2 pairs of purely imaginary eigenvalues, ±icjfc. A; = 3, . . . , n 
Informally, an index two saddle on a potential surface corresponds to a maximum or "hilltop" 
in the potential 
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A variety of algorithms exist for determining critical points on potential surfaces 



33, 



34l |. The number of higher index critical points on a potential surface must be consistent 
with the Morse inequalities (provided the critical points are non-degenerate) 32|, |33| , and 
this constraint is often useful in ensuring that all critical points have been located. 

It has however been argued that critical points of index 2 and higher are of no direct 



chemical significance 



32 



35| . According to the Murrell-Laidler theorem, if two minima on 



a potential surface are connected by an index two saddle, then there must be a lower energy 
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path between them involving only ordinary transition states (index one saddles^ 



(There are certain well-understood limitations to the Murrell-Laidler theorem 37|.) It is 
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therefore true that, provided the reaction coordinate is defined as a minimum energy path 



(MEP) 3^ , it cannot pass through an index two saddle on a potential surface 33| . If however 
;he reaction coordinate is instead associated with a gradient path (steepest descent /ascent) 
391], then it is possible for an index two saddle to lie directly on the reaction path 
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Index two saddles (SP 2) were studied in some detail by Heidrich and Quapp 31i, who 
identifed two (extreme) classes of SP 2 critical points. The first class, "virtual" saddles of 
index two (V-SP 2), are associated with two essentially independent index one transition 
states in different parts of a large molecule, and so can be thought of simply as the direct 
product of two index one saddles; an example of a V-SP 2 saddle occurs in the effectively 
independent internal rotation of the methyl groups in dimethyl ether 4l|. In the second 
class of saddles, "proper" saddles of index two (P-SP 2), the motions associated with the 
two imaginary vibrational frequencies at the saddle ( "downhill" motions) involve the same 
sets of atoms, and so represent inherently two dimensional dynamical processes. 

Examples of proper index two potential saddles abound. Heidrich and Quapp discuss 
the case of face protonated aromatic compounds, where high energy saddle points of index 
two prevent proton transfer across the aromatic ring, so that proton shifts must occur at 
the ring periphery 3l|. Index two saddles are found on potential surfaces located between 
pairs of minima and index one saddles, as in the case of internal rotation/inversion in the 
H2BNH2 molecule 42], or connected to index one saddles connecting four symmetry related 
minima, as for isomerization pathways in B2CH4 43|]. In general, index two saddles are found 



separating symmetry related transition states in so-called "narcissistic" transformations |44]. 

It is quite plausible that, for low enough potential barriers 45|] or at high enough en- 
ergies 461], saddles with index > 1 might well play a significant role in determining sys- 
tem properties and dynamics. There is for example a continuing debate concerning the 
role of higher index saddles in determining the behavior of supercooled liquids and glasses 
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551 . Computations on 256 atom LJ binary mixtures [5 



show that a significant fraction of saddles with index > 2 in hig. 
are however "virtual" saddles in the sense of Heidrich and Quapp 



i-dimensional systems 



3l|. The general rela- 



tion between configuration space topology (distribution of saddles) and phase transitions 
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great current interest 
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56|- Several examples o: 



■ non-MEP or non-IRC reactions 



621 ] and "roaming" mechanisms 
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67|, 



68| have been iden- 



tified in recent years; the dynamics of these reactions is not mediated by a single conventional 
transition state associated with an index one saddle. 

Earnshaw's theorem j69| suggests that higher index saddles are likely to be ubiquitous 



in classical models of atoms and molecules, where all particles interact via the Coulomb 
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potential. The role of index two saddles in the (classical) ionization dynamics of the He 
atom has in fact recently been studied from a phase space perspective by Haller et al. 

Phase space structures and their infiuence on phase space transport were developed in 
some detail for the case of an index one saddle of the potential energy surface (henceforth, 
"index one saddles") corresponding to an equilibrium point of saddle-center-. . .-center sta- 
bility type for the corresponding Hamilton's equations [lO|. In the present paper we follow a 
similar path for index two saddles. (A generalization of phase space transition state theory 
from index one to higher index saddles has independently been given by Haller et al. 7o|].) 

The outline for this paper is as follows. We begin in Sec. [Ill by discussing the phase space 
structure of the inverted parabolic barrier, an index one saddle. This elementary example 
provides the foundation for our treatment of higher index saddles. In Section [XTTl we consider 
geometry and transport associated with a 3 DoF quadratic Hamiltonian for a saddle- saddle- 
center. This simple model contains the essential features of the problem we are considering. 
The quadratic Hamiltonian is separable (and integrable), and we can therefore consider 
the saddle-saddle dynamics and center dynamics separately. The saddle-saddle dynamics 
is considered in detail in Section IIII Bl where we show that the two integrals associated 
with the saddle-saddle dynamics serve to classify the geometry of trajectories that pass 
through a neighborhood of the saddle-saddle. This "reduced dimensionality" analysis of the 
saddle-saddle dynamics associated with separability of the quadratic Hamiltonian is central 
to our analysis. We introduce a symbolic representation of the qualitatively distinct classes 
of trajectory behavior in the vicinity of the saddle-saddle equilibrium, and show that this 
symbolic representation can be placed into direct correspondence with the different classes 
of trajectory in the vicinity of the potential hilltop. 

For completeness, and in analogy with the 3 DoF case, we consider a quadratic n DoF 
model of a Hamiltonian having an equilibrium point of saddle-saddle-center-. . .-center equi- 
librium type in Appendix El and show that separability of the quadratic Hamiltonian implies 
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that the 3 DoF analysis easily goes though for the n DoF case, with obvious modifications. 
In Section HVl we consider the general case of a n DoF fully nonlinear Hamiltonian having an 
equilibrium point of saddle-saddle-center-. . .-center stability type. In a neighborhood of the 
equilibrium point the Poincare-Birkhoff normal form theory can be applied to construct a set 
of new coordinates in which the Hamiltonian assumes a simple form (which we describe). 
In fact, for generic non-resonance conditions on the eigenvalues of the matrix associated 
with the linearization of Hamilton's equations about the equilibrium point, the normal form 
Hamiltonian is integrable. We show that integrability provides all of the advantages that 
separability provided for the quadratic Hamiltonians: the saddle-saddle dynamics can be 
described separately and the two integrals associated with the saddle DoFs can be used to 
characterize completely the geometry of trajectories passing through a neighborhood of the 
saddle-saddle-center-. . .-center equlibrium point. 

As for the case of index one saddles, normally hyperbolic invariant manifolds (NHIMs) 
[sl, [tiI] associated with index two saddles are an important phase space structure and we 
describe the nature of their existence and the role they play in phase space transport in the 
vicinity of index two saddles in Sec. |Vl 

In Section IVl] we discuss two examples where phase space structure in the vicinity of an 
index two saddle is potentially of importance in a problem of chemical dynamics. The first 
example involves isomerization on a potential energy surface with multiple (four) symmetry 
equivalent minima; the dynamics in the vicinity of the index two saddle enables a rigorous 
distinction to be made between stepwise (sequential) and concerted ("hilltop crossing") 
isomerization pathways |45|]. In the second example, two potential minima are connected by 
two distinct transition states associated with conventional index one saddles, and an index 
two saddle sits between the two index one saddles. (The transition states do not have to be 
symmetry equivalent 4^.) For high enough energies, analysis of the phase space geometry 
in the vicinity of the index two saddle enables different classes of reactive trajectories to 
be distinguished. In particular, for the case of non-equivalent index one saddles, we can 
rigorously define "non-MEP" or "roaming" events. Section IVIII concludes with a discussion 
of some outstanding problems for future investigation. 
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II. INDEX ONE SADDLE: INVERTED HARMONIC OSCILLATOR 



In this section we give a quantitative description of trajectories passing through a neigh- 
borhood of an equihbrium point of saddle stabihty type (in brief, a "saddle") for 1 DoF 
Hamiltonian systems. While the basic material is of course well known, our discussion will 
illustrate the distinction between physical coordinates and normal form coordinates, show 
how an integral is used to distinguish qualitatively different trajectories that pass through 
a neighborhood of the saddle, and develop a symbolic description that describes the quali- 
tatively distinct trajectories that pass through a neighborhood of the saddle. Establishing 
these ideas in this simple case will allow for a straightforward generalization to more DoF, 
beginning in Section IIIII 

We consider the problem of a 1-dimensional inverted harmonic oscillator (negative 
parabolic potential), which in suitably scaled physical coordinates {qi,pi) corresponds to 
the Hamiltonian 



with Ai > 0. In general, a Hamiltonian of this form provides an approximate description of 
motion in the vicinity of the top of a (generic, nondegenerate) 1-dimensional barrier. The 
canonical transformation {qi,pi) ^— > {qi,Pi), where 




(1) 




(2b) 



(2a) 



transforms the inverted harmonic oscillator Hamiltonian ([T]) into the normal form 



hi = Xipiqi. 



(3) 



The associated equations of motion in normal form coordinates are 




(4a) 



(4b) 



with solutions 




(5a) 
(5b) 
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The phase space origin (gi = 0,pi = 0) is therefore an equihbrium point of saddle type. 
Phase space portraits for the inverted harmonic oscillator in terms of the physical variables 
{qi,pi) and the normal form variables {qi,Pi) are shown in Figured! We remark that the 
canonical (or symplectic) transformation to normal form coordinates, and the Hamiltonian 
expressed in these coordinates (known as the "normal form Hamiltonian" ) are trivial in this 
1 DoF case. The real power of the normal form approach becomes clear for 2 and more DoF 
where it is a simple matter to describe trajectories near a higher dimensional "saddle" and 
the phase space structures near the saddle can be easily described in terms of (approximate) 



integrals and the normal form coordinates (see 17| and refs therein). 



The Hamiltonian can be written in terms of the action variable 

h = PiQi (6) 

which is obviously a conserved quantity (proportional to the Hamiltonian hi). The action 
Ii can be positive, negative or zero. 

The condition Ji = defines 2 co-dimension one [3] invariant manifolds {pi = and 
qi = 0) that intersect at the saddle equilibrium point, which is the most basic example of 
a NHIM. More importantly, these invariant manifolds divide the phase plane into 4 regions 
(quadrants) corresponding to qualitatively different types of trajectory. 

In terms of the physical variables (gi,pi), the two quadrants with Ji > are associated 
with motions where the particle has enough energy to surmount the top of the parabolic 
barrier, so that there is no turning point in the coordinate qi (classical barrier transmission), 
whereas the two quadrants for which the action Ji < are associated with motions for which 
the particle has insufficient energy to surmount the barrier, so that the coordinate qi exhibits 
a turning point (classical reflection from barrier). 

To obtain a more precise classification of the possible motions in the vicinity of the 
equilibrium point, we define a phase space neighborhood Af of the equilibrium point as 
follows: 

^ = {{qi,Pi) I \qi\ < £1, bil < &i} , (7) 



for some suitably chosen Si [73|. A trajectory enters A/" if at some time it crosses the 
boundary of A/" with the velocity vector pointing into A/". Similarly for trajectories that exit 
Af. It is easy to see from ([5]) that, with the exception of the (zero Lebesgue measure) set 
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of trajectories that enter TV with \qi\ = 0, all trajectories that enter J\f, exit TV. In fact, we 
can characterize the condition for entry and exit of trajectories to and from TV as follows: 

Condition for Entry into TV: |pi(0)| = ei, < e\. 

Condition for Exit from TV: |gi(0)| = £i, |Ji| < e\. 

Trajectories that exit and enter TV remain in the same quadrant in the plane, and 

either pass over or are reflected by the potential barrier (this latter interpretation is based 
on the Hamiltonian expressed in terms of the physical coordinates). Corresponding to the 
four distinct quadrants, we define four classes of trajectory: 

• (++) (91 > 0, pi < 0) 

• (-+) (gi < 0, < 0) 

• (+-) (gi > 0, pi > 0) 

• (— ) (gi < 0, Pi > 0) 

This description uses both the physical and normal form coordinates, and the notation 
indicates the net result of passage through region TV in terms of the action of an associated 
"scattering map". Considering the first item, reading from right to left [tJ], (+, +) refers 
to a trajectory in the physical coordinates that enters the neighborhood of the saddle with 
gi > and exits the neighborhood of the saddle with q\ > 0, i.e. it is reflected from the 
barrier on the qi > side. In the normal form coordinates this trajectory corresponds to a 
trajectory in the fourth quadrant, [qi > 0, pi < 0). Similarly, trajectories of type ( — h) enter 
a neighborhood of the saddle with gi > 0, pass over the saddle, and exit the neighborhood 
of the saddle with qi < 0. In normal form coordinates these correspond to trajectories in 
the third quadrant, (gi < 0, pi < 0). The remaining two cases are understood similarly. 
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III. PHASE SPACE GEOMETRY AND TRANSPORT ASSOCIATED WITH A 
3 DOF QUADRATIC HAMILTONIAN OF A SADDLE-SADDLE-CENTER EQUI- 
LIBRIUM 

In this section we describe the nature of trajectories as they pass through a neighborhood 
of an index 2 saddle for a 3 DoF quadratic Hamiltonian. We will express the Hamiltonian 
in normal form coordinates. However, the symbolic representation of the trajectories will be 
based on dynamics described in terms of the physical coordinates, as in the previous section. 

A. Quadratic saddle-saddle-center system 

We now consider a 3 DoF quadratic Hamiltonian that describes the dynamics near a 
saddle-saddle-center stability type equihbrium point located at the origin: 

Hssc = Ai^ipi + X2q2P2 + ^ (?3 + pT) ' Ai, A2, u; > 0. (8) 

Note that the variables {qi,Pi, q2,P2) in Hamiltonian ([8]) are not physical coordinates, but 
rather normal form coordinates as defined in eq. ([2]). The Hamiltonian ([8]) can be used 
to describe dynamics in the vicinity of a hill-top (nondegenerate index two saddle) in a 3 
dimensional potential energy surface. 

Hamiltonian is completely integrable, with (independent) integrals: 

h = h = q2P2, ^3 = ^ (^3 +P3) > (9) 
and can be expressed as a function of the integrals: 

Hssc = Ai/i + A2/2 + ujh. (10) 

These integrals play a crucial role in our understanding of the underlying phase space ge- 
ometry and transport. 

B. Phase Space Geometry and Transport Associated with the 2 DoF Saddle- 
Saddle Subsystem 

An important consequence of the separability of the 3 DoF quadratic Hamiltonian ([8]) is 
that an understanding of the nature of the trajectories passing through a neighborhood of 
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the saddle-saddle-center equilibrium point can be obtained through a study of the 2 DoF 
subsystem corresponding to the saddle- saddle, since motion in the center DoF is bounded, 
and the action J3 is conserved. The Hamiltonian for this 2 DoF subsystem is: 

Hss = Al^lPl + X2q2P2, (11) 

with integrals Ji and I2 as given in ([9j). 

Hamilton's equations for the 2 DoF saddle-saddle subsystem are given by: 

dpi 

Pi = = -XiPi, (12b) 

oqi 

dHgg . . 

h = — = >^2q2, (12c) 
dp2 

P2 = — ^ = -A2P2 (12d) 
dq2 

The phase flow in the vicinity of the saddle-saddle equilibrium is shown in Fig. [21 In 
this figure we show the phase flow in both the physical coordinates defined in ([1]) and 
the normal form coordinates. The transformation between these two sets of coordinates is 
given in ([2]). Since our Hamiltonian is separable, it is a trivial matter to use this coordinate 
transformation on both saddles to express the Hamiltonian in physical coordinates. However, 
the separability also means that we can immediately adapt the discussion of the previous 
section for physical coordinates to the saddle-saddle case, and it is this approach that we 
will follow. 

As before, we define a neighborhood A/" of the equilibrium point in phase space as follows: 
■Af = {iqi,Pi,q2,P2) I \qi\ < ei, \pi\ < eu \q2\ < 62, \P2\ < £2} , (13) 



for some suitably chosen ei and €2 [73|. As for the 1 DoF case, with the exception of the 
(zero Lebesgue measure) set of trajectories that enter Af with = \q2\ = 0, all trajectories 
that enter A/", exit A/". We characterize the condition for entry and exit of trajectories to 
and from J\f as follows: 

Condition for Entry into A/": |^'i(0)| = £1, |p2(0)| = £2, \hh\ < 6\e\. 
Condition for Exit from M: |q'i(0)| = £1, |q'2(0)| = £2, \^i^2\ < ^fel. 
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As for the 1 DoF unstable equilibrium, we can define different classes of trajectory char- 
acterized by their behavior under the scattering map corresponding to passage through the 
region M . For the 2 DoF case, there are 4x4 = 16 combinations of quadrants, and 16 types 
of trajectory. The symbolic description of the behavior of a trajectory as it passes through 
a neighborhood of the saddle-saddle with respect to the physical coordinates Qk-, k = \,2/is 
expressed by the following four symbols, (/1/2; ^1^2), where ii = ±, 22 = ±, /i = ±, /2 = ±- 
Here ik-, k = 1,2 refer to the "initial" sign of ^a,-, ^ = 1, 2 as it enters the neighborhood 
of the saddle-saddle and fk, k = 1,2 refer to the "final" sign of g^, /c = 1, 2, as it leaves 

the neighborhood of the saddle. For example, trajectories of class ( ; H — ) pass over the 

barrier from gi > to gi < 0, but remain on the side of the barrier with q2 < 0. 

Of the 16 qualitatively distinct classes of trajectory, the 4 types (++; ), ( — h; H — ), 

(H — ; — h) and ( ; -|--|-) are the only trajectories that undergo a change of sign of both 

coordinates qi and q2, and so they are the trajectories that pass "over the hilltop" in the 
vicinity of the saddle-saddle equilibrium. 

As in the 1 DoF case, co-dimension one surfaces separate the different types of trajectory. 
These are the four co-dimension one invariant manifolds given by qi = 0, pi = 0, q2 = 0, P2 = 
(i.e. Ii = 0, I2 = 0). For example, the co-dimension one surface gi = forms the boundary 
between trajectories of type (++; H — ) and ( — h; H — ), and so on. 

The dynamical significance of this symbolic classification is disussed in Section IVT] in the 
context of isomerization reactions. 

C. Including the Additional Center DoF 

Because the system is separable, including the additional center DoF has no effect on our 
discussion above on the nature of trajectories in the 2 DoF saddle-saddle subsystem as they 
pass through a neighborhood of the saddle. 
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To see this, note that equations of motion for Hamiltonian (fTT]) are given by: 

Qi = = Aigi, (14a) 



dp 



= _^^£££ ^ -A^pi, (14b) 
dqi 

h = = ^2^2, (14c) 

p2 = — = -A2P2, (14d) 

dHssc I \ 

q-s = = ^P3, (14e) 

dps 

dHssc r-i AC\ 

P3 = ^ = -W?3, (14f) 

dq3 

As previously, we define a neighborhood of the saddle-saddle-center equilibrium point in 
phase space as follows: 

-A/" = {(gi,Pl,?2,P2,?3,P3) I < £1, \Pl\ < £1, \q2\ < £2, \P2\ < £2 ^31 < £3, b3| < £3} , 

(15) 

for suitably chosen £i, £2 and £3. We want to describe the geometry associated with tra- 
jectories that enter and leave A/". First, note again the fact that Hamilton's equations are 
separable, and that the (simple harmonic) motion of ^3 and p^, is bounded. It then follows 
that if qs and ps are initially chosen to satisfy the condition for being in Af, then they sat- 
isfy this condition for all time. Hence, the issue of trajectories entering Af and exiting Af 
depends entirely on the behavior of the {qi,Pi,q2,P2) components of a trajectory and, as a 
consequence of separability, this behavior is exactly as described in Section IIII B[ 

We now examine some further aspects of the geometry for the 3 DoF quadratic Hamilto- 
nian model of a saddle-saddle-center. The equilibrium point is located at the origin and has 
zero energy. We consider geometric structures in the energy surface for positive energies: 

Hssc = XiqiPi + \2q2P2 + I [ql + P3) = > 0. (16) 

It is clear from f[T^ that = = g2 = P2 = is a two dimensional invariant manifold in 
the 6 dimensional phase space. Using (ITBl) . its intersection with the 5 dimensional energy 
surface is given by: 

^{ql+pf)=E>Q. (17) 

This is a periodic orbit in the energy surface but, more generally, it is an example of a 
NHIM. The coordinates qi,Pi,q2,P2 can be viewed as coordinates for the normal directions 

14 



of the NHIM, and it follows from the nature of the asymptotic (in time) behavior of these 
coordinates (see f|T^ ) that this NHIM has a 3 dimensional stable manifold and a 3 dimen- 
sional unstable manifold in the 5 dimensional energy surface. Further discussion of the role 
of NHIMs is given in Section IV CI 
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IV. n DOF, HIGHER ORDER TERMS IN THE HAMILTONIAN, AND THE 
POINCARE-BIRKHOFF NORMAL FORM 



The examples we have considered so far have been exceptional - quadratic, separable ~ 
Hamiltonians (which implies that they are completely integrable). Now we will show that 
for a general time-independent Hamiltonian, in a neighborhood of an equilibrium point of 
saddle-saddle-center-. . .-center stability type, a "good set of coordinates" can be found in 
terms of which the Hamiltonian assumes a simple form for which the results described above 
hold. Moreover, these coordinates are obtained via a constructive algorithm. 

For completeness, in Appendix|X]we consider in detail the nature of trajectories associated 
with a quadratic n DoF Hamiltonian model that describes the dynamics near a saddle-saddle- 
center-. . .-center stability type equilibrium point, and then we consider the case where higher 
order terms are included in this model. 



A. Normal form for n DoF 

Specifically, in a neighborhood of an equilibrium point of saddle-saddle-center-. . .-center 
stability type, Poincare-Birkhoff normal form theory can be used to (formally) construct a 
symplectic change of coordinates so that in the new coordinates the Hamiltonian has the 
form: 

n 

H = Xiqipi + \2q2P2 + X] y {(fi+Pi) + • • • ,?n,Pl, • • • ,Pn), (18) 

i=3 

where /(<?i, . . . , g„,,pi, . . . ,p„) assumes a particularly simple form that is amenable to anal- 
ysis. This is discussed in some detail, with particular relevance to reaction dynamics, in 



ref. 171]. For our purposes, if we assume that the purely imaginary eigenvalues satisfy the 
non- resonance condition k-^oj-i + . . . + A;„a;„ 7^ for any n — 2 vector of integers {k^, . . . , /c„) 
with not a// the /cj = (that is, (A;3, . . . , kn) G Z'"~^ — {0}) and the real eigenvalues satisfy the 
(independent) non-resonance condition kiXi + k2\2 7^ for any 2 vector of integers (/ci, k2) 
with not all the fcj = 0, i = 1, 2, then /(51, . . . , . . . ,p„) can be represented as an even 

order polynomial in the variables 

h = qiPu h = q2P2, h = — [qk+Pk) ^k = ?,,...,n. (19) 
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In other words, we can express the normal form Hamiltonian in this situation as 



H{h,l2,l3,...,In), (20) 



with the associated Hamilton's equations given by: 



. _dH _ dHdh _ dH 
'^^ dpi dli dpi dli ' 
• _ _9H _ _dH dh _ _dH 

dqi dli dqi dli^^' 
. dH dHdh dH .^^ , 

dH dHdh dH 
dq2 dh dq-2 dh 

. dH dHdh dH . . 

1k = ir- = -WTTT- = TrrPk. (21e) 



dpk dh dpk dh 
dH _ dH dh _ 

dqk dh dqk dh ' 
and it can be verified by a direct calculation that the Ik, k = 1, . . . ,n are integrals of the 



R- = = = -^«?fc, /c = 3, (21f) 



motion for fl2T|) 

Since the h, k = 1, . . . ,n are constant on trajectories, it follows that the partial deriva- 
tives 1^, k = 1, . . . ,n, are also constant on trajectories. Hence, it follows that the structure 
of eqs ( I2T]) is extremely simple, and this is a consequence of the integrability of the normal 
form. 

We can define a neighborhood of the equilibrium point exactly as in eq. flA3l) in Ap- 
pendix |A] and consider the geometry associated with trajectories that enter and leave C. 
Even though the nonlinear Hamilton's equations (12T|) are not separable, it follows from inte- 
grability and the structure of (|2T1) that the motion of qk and pk is bounded, k = 3, . . . ,n { qk 
and Pk just undergo periodic motion, individually for each k, although the frequency of mo- 
tion can vary from trajectory to trajectory). It then follows that if qk and pk, k = 3. . . . ,n, 
are initially chosen in C, then they remain in C for all time. Hence, exactly as in the 
quadratic case, the issue of trajectories entering C and exiting C depends entirely on the 
behavior of the {qi,Pi,q2,P2) components of a trajectory and, as a consequence of integra- 
bility, this behavior is exactly as described in Section IIIIBI and we have the integrals h and 
h as parameters for describing the geometry of the trajectories during their passage though 
C. 
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Exactly as for the model quadratic Hamiltonians, it is simple to verify that qi = pi = 
Q'2 = P2 = is a 2?i — 4 dimensional invariant manifold in the 2n dimensional phase space. 
For energies greater than zero (the energy of the saddle-saddle-center-. . .-center), in the 
2n — 1 dimensional energy surface it is a 2n — 5 dimensional NHIM with 2n — 3 dimensional 
stable and unstable manifolds (they are co-dimension two [t^ in the energy surface). An 
important property of NHIMs is that they, as well as their stable and unstable manifolds, 
persist under perturbation. If we consider a small C, then the terms f{qi,Pi, ■ ■ ■ , QniPn) can 
be considered as a perturbation of the quadratic Hamiltonian, for which the NHIM had the 
geometrical structure of S*^*^"^. If the perturbation is "sufficiently small" then this spherical 
structure is preserved. 



B. Accuracy of the Normal Form 

The normalization proceeds via formal power series manipulations whose input is a Tay- 
lor expansion of the original Hamiltonian, if, necessarily up to some finite order, M, in 
homogeneous polynomials. For a particular application, this procedure naturally necessi- 
tates a suitable choice of the order, M, for the normalization, after which one must make a 
restriction to some local region, £, about the equilibrium point in which the resulting com- 
putations achieve some desired accuracy. Hence, the accuracy of the normal form as a power 
series expansion truncated at order M in a neighborhood C is determined by comparing the 
dynamics associated with the normal form with the dynamics of the original system. There 
are several independent tests that can be carried out to verify accuracy o 



Straightforward tests that we have used are the following 
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14 
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the normal form. 
3, El: 



• Examine how the integrals associated with the normal form change on trajectories of 
the full Hamiltonian (the integrals will be constant on trajectories of the normal form). 

• Check invariance of the different invariant manifolds (i.e. the NHIM and its stable 
and unstable manifolds) with respect to trajectories of the full Hamiltonian. 

Both of these tests require us to use the transformations between the original coordinates 
and the normal form coordinates. Specific examples where M, C and accuracy of the normal 
forms and the constancy of integrals of the truncated normal form are considered can be 
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found in 111 
found in ll 



13 



14 



15 



771]. A general discussion of accuracy of the normal form can be 
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V. COMPARISON OF PHASE SPACE GEOMETRY AND TRAJECTORIES FOR 
INDEX ONE AND INDEX TWO SADDLES 

In this section we discuss a number of issues that arise from the above results and dis- 
cussions. For comparative purposes it is helpful first to summarize the nature of trajectories 
that pass near an index one saddle and trajectories that pass near an index two saddle. 



A. Index one saddles 



For index one saddles it was shown in previous work that for energies above that of 



the saddle-center-. . .-center there exists a 2n — 3 dimensional normal 
manifold (NHIM) in the 2n — 1 dimensional energy surface j^, 10 . 



y hyperbolic invariant 



7l|. The NHIM is of 



saddle stability type and has 2n — 2 dimensional stable and unstable manifolds. Since these 
are one less dimension than the energy surface, they act as higher dimensional separatrices 
and serve to partition phase space into regions corresponding to qualitatively different types 
of trajectories: 

• Forward reactive trajectories = (H — ) 

• Forward non-reactive trajectories = ( ) 

• Backward reactive trajectories = ( — h) 

• Backward non-reactive trajectories = (++) 

where we have indicated the correspondence with the trajectory classification introduced in 
SecH 

The NHIM is the "anchor" for the construction of a phase space dividing surface having 
the geometrical structure of a 2?i — 2 dimensional sphere. The NHIM is the equator of 
this 2n — 2 dimensional sphere, hence dividing it into two hemispheres. Forward reacting 
trajectories pass through one hemisphere and backward reacting trajectories pass through 
the other hemisphere. These "dividing hemispheres" have the no-recrossing property in 
the sense that the vector field defined by Hamilton's equations is strictly transverse to each 
hemisphere. Moreover, the magnitude of the flux is minimal (and equal) for each hemisphere. 

The language used for the description of trajectories for the case of an index one saddle 
is that of reaction. Trajectories evolve from reactants to products by passing through a 
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dividing surface that locally "divides" the energy surface into two disjoint pieces-reactants 
and products. 

Here, we have classified trajectories in the vicinity of an index one saddle in a slightly 
less general but more explicit fashion in terms of their behavior with respect to crossing a 
potential barrier; this classification is useful in the treatment of reactive dynamics in the 
vicinity of higher index saddles on potential energy surfaces, as discussed below. 

B. Index two Saddles 

The geometrical structures associated with index one saddles provide a rigorous way 
of partitioning the phase space into regions that describe the evolution of a system from 
reactants to products by passing through a dividing surface having the no-recrossing and 
minimal flux properties. This partitioning is made possible by the existence of invariant 
manifolds that have the proper dimension and stability type. Moreover, their geometric 
character is such that the partitioning that they provide has the natural interpretation 
of regions of reactants and products of the energy surface, and the stable and unstable 
manifolds of the NHIM provide natural boundaries between these regions. 

For index two saddles the same types of invariant manifolds still exist - a NHIM and co- 
dimension one (in the energy surface) invariant manifolds. These co-dimension one invariant 
manifolds (g^ = 0, pi = 0, i = 1,2) were introduced above and are discussed further below. 
It was shown that the passage of trajectories through the neighborhood of the equilibrium 
could be understood in terms of a 2 DoF saddle-saddle subsystem and values of the two 
associated action integrals. 

For index one saddles not all trajectories could pass through the dividing surface (by 
definition, the forward and backward non- reactive trajectories do not do so). However, for 
index two saddles all trajectories (except for the set of Lebesgue measure zero discussed 
above) that enter the neighborhood of the equilibrium point exit the neighborhood of the 
equilibrium point, and there are sixteen qualitatively different types of trajectory that enter 
and exit this neighborhood. These sixteen classes are characterized by their crossing behavior 
with respect to the "hilltop" in the potential surface associated with the index two saddle. 

Two examples involving isomerization reactions where this classification of trajectories is 
of chemical relevance are discussed in Section IVIi 
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C. The Role of NHIMs and the existence of co-dimension one invariant manifolds 
in the energy surface 



It has been previously been argued 10|] that the invariant manifold structure associated 
with index one saddles is a natural generalization to 3 or more DoF of the PODS, and their 
associated invariant manifolds, where the latter have provided a fundamental understanding 
of transition state theory for systems with 2 DoF |3|. The necessary ingredient to make 
the leap from 2 DoF to higher DoF is the realization that a NHIM of saddle stability type 
is the natural generalization to higher DoF of the periodic orbit of saddle stability type 

The resulting geometric picture of reaction for N > 2 DoF is 



familiar from 2 DoF 



10 



essentially that of the 2 DoF case, the key feature of the generalization being the availablity 
of the necessary invariant manifolds to partition the energy surface in the multidimensional 
case in a manner that is meaningful for reaction. 

The situation with index two saddles is fundamentally different. On energy surfaces 
above that of the saddle, a NHIM still exists, but it is co-dimension 4 in the energy surface 
with co-dimension 2 stable and unstable manifolds. Nevertheless, as shown above, it is 
still possible to construct co-dimension one invariant manifolds (in the energy surface) that 
locally partition the energy surface into regions associated with qualitatively different types 
of trajectory passing through a neighborhood of the equilibrium point: the surfaces defined 
by the conditions qi = 0, pi = q2 = 0, and P2 = are each co-dimension one invariant 
manifolds in the 2n dimensional phase space, and they are also co-dimension one invariant 
manifolds when intersected with the energy surface. They do not have the property of being 
compact and boundaryless like the NHIMs we have discussed above. We would not expect 
this since (almost) all trajectories that enter C leave C These co-dimension one invariant 
manifolds intersect the boundary of C, and their projection into the qi—pi—q2—P2 space gives 
the coordinate axes shown in Fig[2l This is consistent with our description of trajectories 
that enter and exit C through the use of the integrals Ii and l2- 

While the existence of co-dimension one invariant manifolds in phase space associated 
with index one or index two saddles is conceptually appealing, practically, their main use 
derives from the integrable structure of the normal form and the expression of Hamilton's 
equation using this integrability (e.g. (!2T]) ). This expression allows for a decoupling of 
the dynamics into unbounded motion (the "reactive modes" described by /i and I2) and 
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bounded motion (the "bath modes" described hj Ik, k = 3, . . . ,n). The "reactive dynamics" 
is described by a reduced dimensional system and the dynamics in this system is integrable 
and characterized by the level sets of these integrals. The relevant co-dimension one invariant 
manifolds drop out of this description naturally. 

In the earlier work related to index one saddles in j^, [lo| the integrable structure was 
not emphasized and the uses of the integrals was not developed. Rather, the emphasis was 
on the use of the normal form in discovering the phase space structures. It was hoped that 
a by-product of the normal form computations might be useful and intuitive visualization 
methods of the kind that have proved extremely useful and insightful for 2 DoF systems. 
However, visualization schemes for the invariant manifolds gover ning reaction in systems 



with three and more DoF have thus far been of limited value. In [13|] attempts were made 
to visualize the phase space invariant manifolds governing (planar) HCN isomerization in 
three dimensions by projecting the manifolds into the three dimensional configuration space. 
In this case, the dividing surface (in a fixed energy surface) was four dimensional with its 
equator, the NHIM, three dimensional. It was difficult to interpret their projections into the 
three dimensional configuration space (visually, they did not appear significantly different). 
Moreover, reversibility of Hamilton's equations in this case implies that the projections of 
the four dimensional stable and unstable manifolds into configuration space are identical. 
While the visualization of high dimensional phase space structures is an appealing concept, 
we believe that the dimensional reduction afforded by integrability of the normal form in 
a neighborhood of the equilibrium point is a much more powerful tool for understanding 
reactive phenomena. It provides a precise notion of dimensional reduction to the number 
of DoFs undergoing "reactive behaviour" (i.e. the DoF for which there is the possibility 
of leaving a neighborhood of the equilibrium point). The significance of integrability was 
emphasized in 17|| where it was also shown how the phase space structures can be expressed 
in terms of the integrals (which is what we have essentially done here). While their utility 
for higher dimensional visualization remains to be demonstrated, for index one saddles the 
divid,.g ..face play. a. i,„po.ta,. .o.e fo. a ce.ai. type, of sampling flU Q and the 
NHIM plays a role in determining the fiux through the dividing surface. We will discuss this 
in more detail in Section IVIIi 
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D. Higher index Saddles 

Once the generalization from index one to index two saddles has been made, generalization 
to higher index saddles is straightforward. For an n DoF deterministic, time-independent 
Hamiltonian system, an index k saddle (1 < A; < n) of the potential energy surface corre- 
sponds to an equilibrium point of Hamilton's equations where the matrix associated with the 
linearization of Hamilton's equations about the equilibrium point has k pairs of real eigen- 
values and n — k pure imaginary pairs of eigenvalues. Assuming that appropriate generic 
non-resonance assumptions on the real and imaginary eigenvalues hold, then a transforma- 
tion to an integrable Poincare-Birkhoff normal form in a neighborhood of the equilibrium 
point as described above can be carried out. For energies above the equilibrium point, the 
system has k integrals describing "reactive" dynamics (i.e., trajectories that can enter and 
leave a neighborhood of the equilibrium point) and n — k integrals describing "bath modes" 
(i.e., bounded motions in a neighborhood of the equilibrium point). Hence integrability 
allows the reactive dynamics to be studied with a 2k dimensional reduced system with the 
k reactive integrals providing precise information about the nature of trajectories that pass 
through a neighborhood of the equilibrium point. Moreover, for energies above the saddle, 
the system has a NHIM of dimension 2n — 2A; — 1 having 2n — k — 1 dimensional stable and 
unstable unstable manifolds in the 2n — 1 dimensional energy surface. In the neighborhood 
of the equilibrium point, there are 4}^ qualitatively distinct trajectory classes. 



E. Quantitative Application of Poincare-Birkhoff" Normal Form Theory 

Applying the formalism for higher index saddles to realistic molecular systems requires 
software that computes the normal form Hamiltonian, extracts the integrals of the mo- 
tion, and computes the coordinate transformation (and its inverse) from the original phys- 
ical coordinates to the normal form coordinates, all to sufficiently high order to yield the 
desired accuracy (and accuracy must be assessed by a battery of tests at all points of 
the calculation). Software to carry out this program for index one saddles is available at 



http://lacms.maths.bris.ac.uk/publications/software/index.html In principle it is possible 



to modify this software to carry out the analysis described above for index two, and higher, 
saddles. 
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VI. INDEX TWO SADDLES AND ISOMERIZATION DYNAMICS 



A number of possible contexts in which index two saddles might be important have 
already been discussed in Sec. HI It is however useful to specify the possible dynamical 
role of index two saddles more precisely. We consider two examples: first, the dynamical 
definition of stepwise (sequential) versus concerted isomerization mechanisms on a multi- 
well potential surface; second, the classification of reactive trajectories in a system in which 
there are two distinct transition states. In the latter case our analysis of dynamics in the 
vicinity of the index two saddle suggests a rigorous identification of so-called "non-MEP" 
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62( 1 or "roaming" trajectories 
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68( 1 based on phase space 



structures. 

It should be noted that the partitioning of phase space we discuss is local, that is, re- 
stricted to the vicinity of the index two saddle. 

A. Stepwise vs concerted isomerization 

Index two saddles are of dynamical significance in the identification of sequential versus 
concerted isomerization mechanisms in the case of multiple well isomerizations. 

1. Stepwise versus concerted isomerization in 2 DoF 



Consider a system with two coordinates {qi,q2), having a multiwell potential surface 
v{qi,q2) of the kind shown schematically in Figure [3] [3]. Such a surface can describe. 



for example, conformational 



energies as a function of internal a ng. 
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es 



) 



411 ] or internal rota- 



tion/inversion potentials 42[. (See, for example. Figure 1 in ref. 

The model potential shown has 4 minima. Each potential minimum is given a symbolic 
label according to the quadrant of the (^1,^2) plane in which it is located. For example, 
minimum ( ) is located in the quadrant (gi < 0,^2 < 0). 

The minima are connected by index one saddles, each denoted by the symbol |. At 
higher energy, in the middle of the potential surface, sits an index two saddle, denoted |J 
(the "hilltop"). (We locate the index two saddle at the coordinate origin, (^1,^2) = (0,0), 
and take the value of the potential at the hilltop to be zero, v{0, 0) = 0.) 
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For total energies E < 0, the system can only make transitions between minima by 
passing through the phase space dividing surfaces (transition states) associated with the 
index one saddles (NHIMs). For such energies, we can define the distinct bound (reactant) 
regions of phase space in the usual way, and describe the isomerization kinetics in terms of 
the associated phase space volumes and fiux through phase space dividing surfaces associated 
with the index one saddles. 

Classes of isomerizing trajectories associated with particular transitions between wells 

can be denoted symbolically: for example, trajectories of type ( — h; ) start in the well 

( ) and pass to the well ( — h) through the appropriate dividing surface. 

It is clear that, for E < 0, the only way the system can pass from well ( ) to well (++), 

say, is by a stepwise mechanism: a sequence of isomerizations such as (H — ; ) followed 

by (++; H — ) is required. 

By contrast, for E > 0, trajectories are in principle able to "roam" around on the potential 
surface. In addition to the sequential mechanism just described, there is another possibility, 
namely, a hilltop crossing or concerted mechanism, where, for the example just discussed, 

a trajectory simply passes "directly" from well ( ) to well (++) without undergoing the 

stepwise isomerization process which is necessary for E < 0. 

It is then natural to ask: how can we distinguish between these mechanisms in a dy- 
namically rigorous way? It is not sufficient to analyze trajectories in configuration space, 
as both sequentially isomerizing and hilltop crossing trajectories will pass close to the index 
two saddle. 

The key to dynamical analysis of isomerizing trajectories is the symbolic trajectory clas- 
sification introduced previously. We imagine that the saddle-saddle normal form has been 
obtained in the vicinity of the index two saddle. This means, in particular, that a (nonlin- 
ear, in general) invertible symplectic transformation T has been obtained relating physical 



coordinates {qi,Pi, q2,P2) and normal form coordinates {qi,Pi,q2,P2) [10 



iq,p)=Tiq,p). (22) 

In a phase space neighborhood of the saddle-saddle equilibrium, it is therefore possible to use 
the inverse transformation {q,p) = T~^{q,p) to determine the normal form coordinates along 
a given trajectory. Given these normal form coordinates and the approximate integrability 
of the motion in the vicinity of the index two saddle, a trajectory with energy E > can be 
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classified according to the symbolic scheme developed previously. 

For the particular kind of potential surface considered here, the symbolic classification 
scheme of trajectories using the normal form is precisely the desired classification scheme 

for isomerization behavior. For example, a trajectory of type (H — ; ) passes from well 

( ) to well (h — ), and so is part of a sequential isomerization process, while a trajectory of 

type (++; ) is a hilltop crossing trajectory that, as determined by a rigorous phase space 

criterion in the vicinity of the saddle-saddle equilibrium, passes directly from well ( — ) to 
well (++). 

2. Stepwise versus concerted isomerization in 3 DoF 

We now introduce a third mode: for simplicity, we first consider an uncoupled harmonic 
oscillator DoF. As discussed in Sec. IIII C\ there is now a saddle-saddle-center equilibrium. 
For given total energy E, there will be distinct regimes, depending on the partitioning of 
energy between mode 3 and the 2 DoF isomerizing subsystem. 

If the total energy E is above that of the index two saddle, but there is sufficient energy 
in mode 3, then hilltop crossing isomerization is not possible in the 2 DoF subsystem. If 
on the other hand the amount of energy in mode 3 is small, then the 2 DoF subsystem 
can undergo both sequential and hilltop crossing isomerization. Isomerizing trajectories in 
the 3 DoF system can then be classified using the normal form including the center DoF as 
developed in Sec. IIII CI 

Suppose that mode 3 and the 2 DoF subsystem are now coupled. If we consider the 
behavior of trajectories initialized in a region of phase space with a large amount of energy 
in mode 3 a relevant question is: what is the mechanism by which the system will pass from 
one potential minimum to another, given that the minima in question are separated by more 
than one saddle at low energies? 

Intuitively, there are two possibilities: if significant energy transfer between mode 3 and 
the 2 DoF subsystem does not occur, then the system simply passes from one well to the 
other in a sequential fashion. If however energy transfer occurs between mode 3 and the 
2 DoF isomerizing subsystem, it is possible that sufficient energy can pass into the 2 DoF 
subsystem so that trajectories can isomerize via the hilltop crossing mechanism. In this 
case, hilltop isomerization is mediated by the additional mode. 
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These dynamical alternatives can in principle be distinguished using the transformation 
to normal form coordinates in the vicinity of the saddle-saddle-center equilibrium (cf. Sec. 

Emi). 



B. Competing transition states and non-MEP mechanisms 



We now discuss a second example where analysis of the phase space structure in the 
vicinity of an index two saddle is potentially useful for obtaining deeper mechanistic under- 
standing of an isomerization reaction. 

Consider the model 2 DoF potential shown in Figure HI This potential has two minima, 
two index one saddles (having different energies) and one index two saddle (energy E = 0). 
The index two saddle sits between the two index one saddles (as per Murrell-Laidler [361). 
(See, for example, the double proton transfer reaction in napthazarin. Scheme 14 in ref. l38l .) 

At total energies E < 0, the system can pass from one well to another through either of 
two transition states; in the example shown, the transition state associated with the index 
one saddle at g2 < has lower energy, and hence is associated with the "minimum energy 
path" mechanism. 

At energies above that of the higher energy TS (denoted channel 2), there may exist some 
kind of dynamical competition between transition states: the lower energy TS (channel 1) 
might be narrow (small flux), while the higher energy TS (channel 2) might be broader 
(larger flux). 

For energies above the hilltop, E > 0, there are in principle three possibilities: trajectories 
can isomerize via channel 1, via channel 2, or they can isomerize by "passsage over the 
hilltop". 

Once again, the symbolic classification of trajectories in the vicinity of the saddle-saddle 
equilibrium enables the three classes of trajectory to be rigorously identified. Trajectories of 

class (H — ; ) correspond to the first mechanism (channel 1), class (++; — +) the second 

mechanism (channel 2), while trajectories of class (H — ; — h) or (++; ) are associated 

with the hilltop crossing isomerization. 

All isomerizing trajectories except those of class (H — ; ) can therefore be labelled "non- 



MEP" trajectories [57 



stories except 

58, 3 SJei 



as "roaming" trajectories 63 



64 



3 



62 



66 



. All such non-MEP trajectories can be regarded 



671, 
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VII. SUMMARY AND CONCLUSIONS 



In this paper we have analyzed the phase space structure in the vicinity of an equihb- 
rium point associated with an index two saddle or "hilltop" on a potential energy surface 
("index two saddle"). For the case of model quadratic Hamiltonians we have shown that the 
behavior of trajectories passing through a phase space neighborhood of the equilibrium is 
fully characterized by the value of associated integrals of the motion. We have introduced a 
symbolic classification of trajectories in the vicinity of the saddle-saddle equilibrium, where 
the symbolic representation of trajectory classes is directly related to the nature of the as- 
sociated hilltop crossing dynamics. For the general case where nonlinear terms are present 
in the Hamiltonian, Poincare-Birkhoff normal form theory provides a constructive proce- 
dure for obtaining an integrable approximation to the full Hamiltonian in the vicinity of the 
equilibrium, allowing a corresponding classification of trajectories to be made in the general 
case. 

As for the case of index one saddles, normally hyperbolic invariant manifolds (NHIMs) 
associated with index two saddles are an important phase space structure, and we have 
described the role they play in phase space transport in the vicinity of index two saddles. In 
particular, we have shown that the normal form transformation enables co-dimension one 
surfaces to be defined in phase space; these surfaces serve to partition phase space in the 
vicinity of the saddle-saddle equilibrium into regions corresponding to distinct trajectory 
types. 

Two examples of the importance of index two saddles in problems of chemical dynamics 
were considered. First, isomerization on a potential energy surface with multiple symmetry 
equivalent minima; second, isomerization in a system having two potential minima connected 
by two distinct transition states associated with conventional index one saddles, with an 
index two saddle situated between the two index one saddles. We suggest that classification 
of different classes of reactive trajectories in the vicinity of the index two saddle enables 
a rigorous definition of non-MEP trajectories to be given in terms of local phase space 
structure. 

To conclude, we list several topics that we believe merit further investigation for the case 
of index two saddles. 

Reaction Coordinates. For the case of index one saddles a natural definition of a phase 
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space reaction coordinate is obtained by putting all of the energy into the reactive 



mode 
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171 ]. i.e., for energy H = E, and the Hamiltonian expressed as a function of 
the integrals, we set 

H{h,0,...,0) = E. (23) 

This condition defines a path in the "reactive plane" defined by qi and pi in normal 
form coordinates, which can then be mapped back into the original physical coordi- 
nates by the normal form transformation. 

For the case of an index two saddle we have two reactive modes, described by the 
integrals Ji and I2, and a four dimensional "reactive space" defined by the coordinates 
Qi — Pi — Q2 — P2- Can dynamically significant paths analogous to that for the index 
one saddle be defined in this reactive space? 

Activated Complex. For index one saddles the NHIM has the interpretation as the "ac- 



tivated complex" id, ll7| . Do the NHIMs associated with index k saddles, 1 < k < n 



have a similar chemical significance? 

Flux. For the case of index one saddles the NHIM is the "boundary" of the forward and 
backward dividing surfaces. The magnitudes of the fiuxes across each of these dividing 
surfaces are equal and, using Stokes theorem, are expressed as the integral of an 
appropriate quantity over the NHIM, equal to the volume in the space of "bath modes" 
enclosed by the contour defined by H{0, J2, ...,/„) = i? [l^ . 

For index two (and higher) saddles it is natural to ask whether or not there is any 
chemical significance to the volume in the space of "bath modes" enclosed by the 
contour defined by H{0, 0, /s, . . . , J„) = E. 



Gap times and sampling. Can one develop a useful gap time formalism [80|, |81|, |82| to 

characterize rates of passage in and out of phase space regions defined by the co- 
dimension one surfaces associated with index two saddles (and higher)? 

Related to the gap time formalism, it has been shown that the dividing surface provided 
by the normal form provides a natural surface on which to sample initial conditions 
for trajectory calculations that yield quanities of physical interest such as the reactant 



volume or dissociation rate 
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82l |. Does the normal form associated with index 



30 



two saddles yield manifolds on which a similar trajectory sampling strategy can be 
implemented? 

The Manifolds gi = 0, = 52 = 0, and p2 = , a new type of NHIM? We have 
been careful not to refer to the co-dimension one invariant manifolds defined by 
qi = 0, pi = q2 = 0, and p2 = as "NHIMs". There are two reasons for this. The 
first is that they are not of the character of the NHIMs associated with index one 
saddles, i.e., compact and boundaryless (diffeomorphic to spheres), having (invariant) 
stable and unstable manifolds (locally) partitioning the energy surface into reactants 
and products, and "persistent under perturbations" (practically, this means that the 
manifolds maintain their character for a range of energy above that of the saddle). 
The second is that they do not immediately satisfy the hypotheses of the existing 



normally hyperbolic invariant manifold theorems 7l|. This latter point deserves more 
discussion. Roughly speaking, the property of "normal hyperbolicity" means that, 
under linearized dynamics, the magnitude of the growth and decay of vectors normal 
to the manifolds dominate the magnitude of the growth and decay of vectors tangent 
to the manifold. The manifolds under consideration are obtained by setting one 
coordinate to zero. Therefore, the growth of vectors under the linearized dynamics 
in the direction of the coordinate set to zero describes the normal dynamics, and it 
is a simple matter to verify that in this direction vectors are either exponentially 
growing or decaying (depending on the manifold under consideration). The problem 
is that there will also be directions tangent to the manifold that are exponentially 
growing or decaying. Additionally, these manifolds are not compact and boundaryless. 
Sometimes this difficulty can be dealt with through a consideration of "overflowing" or 



"inflowing" invariance (see 7ll] for examples and related references). Nevertheless, the 
manifolds do exist for the normal form truncated at any order, and this is sufficient 
for our purposes. Using this property, and the ideas and setting just described, it may 
be possible to prove something like a "persistence theorem for normally hyperbolic 
invariant manifolds" for this setting, but we do not pursue that problem here. 

Addressing several of these issues will require globalizing the invariant manifolds obtained 
in the vicinity of the saddle-saddle equilibribrium. Indeed, globalization constitutes an 
outstanding challenge for the theory of the dynamics associated with higher index saddles. 
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APPENDIX A: PHASE SPACE GEOMETRY AND TRANSPORT ASSOCI- 
ATED WITH AN n DOF QUADRATIC HAMILTONIAN OF A SADDLE-SADDLE- 
CENTER-. . -CENTER EQUILIBRIUM POINT 

The following quadratic Hamiltonian describes the dynamics near a saddle-saddle-center- 
. . .-center stability type equilibrium point located at the origin: 

n 

Hn-quad = ^mVl + ^292^2 + {^i + Pi) > (^1) 



2 

j=3 



with associated Hamilton's equations: 

dHYi—quad 



qi = r^""" = Xiqi, (A2a) 
dpi 

Pi = - ^^r"""' = -AiPi, (A2b) 
oqi 

q2 = ^^r^""' = A2g2, (A2c) 
dp2 

P2 = -^^^ = -A2P2, (A2d) 
dq2 

dH-n-quad ■ o / A o \ 

Qi = —-i — = LOiPi, z = 3, . . . , n (A2e) 

dpi 

dHn-quad • o /A of N 

Pi = 7:^ = -^i^j, « = 3, . . . , n. (A2t) 

dqi 

It is straightforward to verify that Ii = qiPi, I2 = q2P2, Ik = ^ + Pk) , k = 3, . . . ,n are 
(independent) integrals of the motion for (]A2p . As previously, we define a neighborhood of 
the saddle-saddle-center equilibrium point in phase space as follows: 

^ = {iqi,Pi,---,qn,Pn) I \qi\<£i,\pi\<£i,i = '^,---,n}, (A3) 

for appropriately chosen £j > 0, z = 1, . . . , n. Again, the fact that Hamilton's equations are 
separable, and that the motion of qk and pk is bounded, k = 3, . . . ,n means that the issue of 
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trajectories entering L and exiting £ depends entirely on the behavior of the ^'2,^2) 
components of a trajectory and, as a consequence of separabihty, this behavior is exactly as 
described in Section [11161 

As for the 3 DoF case considered in Sec. IIIIl the equilibrium point is located at the origin 
and has zero energy, and we will consider geometric structures in the energy surface for 
positive energies: 

n 

Hn^guad = XiQiPi + \2q2P2 + y [qI + pI) = E > 0. (A4) 

k=3 

It is clear from (]A2p that gi = pi = g2 = P2 = is a two dimensional invariant manifold in 
the 2n dimensional phase space. Using flA4p . its intersection with the 2n — 1 dimensional 
energy surface is given by: 

n 

Y.^{<& + pI)=E>^. (A5) 

fc=3 

This is the equation for a 2n — 5 dimensional sphere, S*^*^"^, in the energy surface; as above, it 
is an example of a NHIM [sl, [tiI. The coordinates qiiPi-, q2,P2 can be viewed as coordinates 
for the normal directions of the NHIM, and it follows from the nature of the asymptotic (in 
time) behavior of these coordinates (see ( 1A2I) ) that this NHIM has a 2n — 3 dimensional 
stable manifold and a 2n — 3 dimensional unstable manifold in the 2n—l dimensional energy 
surface. 
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FIG. 1: The phase space geometry associated with the inverted parabohc potential in (a) 
physical coordinates and (b) normal form coordinates. 
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FIG. 2: The geometry associated with the 2 DoF saddle-saddle subsystem. 
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FIG. 3: Contour plot of a model 2 DoF potential having 4 minima, 4 index one saddles (J) 
and one index two saddle (W). 
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